Optimum pinning of the vortex lattice in extremely type-II layered superconductors 
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The two-dimensional (2D) vortex lattice in the extreme type-II limit is studied by Monte Carlo 
simulation of the corresponding 2D Coulomb gas, with identical pins placed at sites coinciding with 
the zero-temperature triangular vortex lattice. At weak pinning we find evidence for 2D melting into 
an intermediate hexatic phase. The strong pinning regime shows a Kosterlitz-Thouless transition, 
driven by interstitial vortex/anti-vortex excitations. A stack of such identical layers with a weak 
Josephson coupling models a layered superconductor with a triangular arrangement of columnar 
pins at the matching field. A partial duality analysis finds that layer decoupling of the flux-line 
lattice does not occur at weak pinning for temperatures below 2D melting. 
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I. INTRODUCTION 

It well known that the motion of vortex lines in the 
mixed phase of a type-II superconductor generates dissi- 
pation, and hence that an unpinned vortex-lattice state 
is in fact resistive*. This has been confirmed recently in 
the mixed phase of clean high-temperature superconduc- 
tors, where the superconductivity of samples with a strip 
geometry is found to be very much superior to that in 
samples with a Corbino disk geometry 2 ". Surface barriers 
in the strip geometry prohibit rigid motion of the vortex 
lattice, while the Corbino disk geometry allows for rigid 
rotations of the vortex lattice. The above phenomena can 
be understood theoretically in the extreme type-II limit, 
where magnetic screening is absent*. Rigid translations 
of the vortex lattice result in an infrared divergence that 
destroys phase coherence at any temperature*. This in- 
frared divergence can be removed by excluding rigid mo- 
tion of the vortex lattice through surface barriers, where- 
upon phase coherence is restored*. 

Defects in the bulk of a superconductor can also effec- 
tively prohibit the rigid motion of the vortex lattice*. In 
this work we study the nature of phase coherence in an 
extremely type-II layered superconductor, with magnetic 
field oriented perpendicular to the layers, and containing 
an array of correlated pins. We choose to use an opti- 
mum arrangement of identical columnar pins 5 , the loca- 
tions of which match the triangular vortex lattice at zero 
temperature. Such a configuration can be realized exper- 
imentally by artificial "anti-dot" arrays*. The Josephson 
coupling between layers is turned off initially, thus al- 
lowing us to model the system of vortices in each layer 
by a two-dimensional (2D) Coulomb gas with a uniform 
charge background and a commensurate pinning poten- 
tial, (cf. ref. 7) We employ Monte Carlo (MC) sim- 
ulations to uncover the thermodynamic phase diagram 
of this system under periodic boundary conditions. The 
depth of the pinning potential, U, becomes a useful con- 
trol parameter. As U increases from zero, the "float- 
ing" vortex lattice phased that exists at U = becomes 
pinned at a critical U p . A finite-size analysis shows that 
Up tends to zero in the thermodynamic limit, and that 



floating is prohibited initially by a sparse distribution of 
pinned vortices. Phase coherence is then restored at a yet 
stronger pinning, U m . The vortex lattice inside the range 
Up < U < U m therefore shows no phase coherence de- 
spite the fact that it does not float! The identification of 
this regime with the hexatic phase 10 is indicated by recent 
theoretical work*. Indeed, a modest size analysis demon- 
strates that this phase shows (strict) long-range orienta- 
tional order, which is a hallmark of the hexatic phase* -. 
The strong-pinning regime at the other extreme exhibits 
a standard Kosterlitz-Thouless (KT) transition driven by 
the unbinding of vortex/anti-vortex pairs that are not 
linked to the vortex latticei*, which is now fixed to the 
pins and appears to be irrelevant** 2 ^ Lastly, the effects 
of a weak Josephson coupling between layers is deter- 
mined through the application of a partial duality analy- 
sis of the corresponding layered XY model 12 . On this ba- 
sis, we conclude that the superconducting-normal tran- 
sition shown by such an optimally pinned vortex lattice 
lies within the universality class of the three-dimensional 
(3D) XY model. 



II. 2D VORTEX LATTICE WITH 
COMMENSURATE PINS 

Consider an infinite stack of isolated superconducting 
layers in a perpendicular external magnetic field. Each 
layer is assumed to be identical in order to reflect the 
correlated pinning. A weak Josephson coupling will be 
switched on later. Magnetic screening effects can be ne- 
glected in the extreme type-II limit assumed through- 
out, in which case the XY model over the square lattice 
with uniform frustration provides a qualitatively correct 
description of the mixed phase of each layer. The cor- 
responding Boltzmann distribution is set by the energy 
functional 
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for the superfiuid kinetic energy in terms of the super- 
conducting phase (j)(r). Here A M 0(r) = 4>(f + aft) — 4>{r) 
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and A — (0,2n fx/a) make up the local supercurrent, 
where / denotes the concentration of vortices over the 
square lattice, with lattice constant a. The local phase 
rigidity J M (r) is assumed to be constant over most of the 
nearest- neighbor links (r, r + a/t), with the exception of 
those links in the vicinity of a pinning site. We shall next 



take the Villain approximation which is generally valid at 
low temperatur e) 13 ! 14 . After making a series of standard 
manipulations, we obtain a Coulomb gas ensemble with 
pins that describes the vortex degrees of freedom on the 
dual square lattice. The ensemble is weighted by the 
Boltzmann distribution set by the energy functional 
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in terms of the integer vorticity field Q(R) over the sites 
R of the dual lattice. The logarithmic interaction be- 
tween the vortices can be expressed as a Greens func- 
tion JG^ = ^2 n \n)e n (n\, where the states \n) diag- 
onalize the operator —^[J^ 1 ^ + Jy 1 ^] with corre- 
sponding eigenvalues e^ 1 . Here, the lattice difference op- 
erator Ay( x ) ac ts between those adjacent points on the 
dual lattice that are split by the link on the base lat- 
tice that the coupling constant J x (y) refers to. The ef- 
fective coupling energy J is set by the requirement that 
(qi\G^ I92} = ?j~ a~ 2 #i,2 for plane waves \<fi) in the long 
wavelength limit, qi — > 0. The pinning potential, on the 
other hand, originates from the contribution to the for- 
mer Greens function by bound states localized at the pin- 
ning sites: V P (R) = (2tt) 2 £ £ „<o e n \{R\n)\ 2 . Finally it is 
instructive to point out that the phase rigidity can also 
be directly computed within the Villain approximation, 
where it is given byi^ one over the dielectric constant of 
the Coulomb gas ensemble {2J. In particular, for anixL 
square mesh with periodic boundary conditions, we have 

p./ J = l-lim(27r) 2 /3 ((%Q_£) - (Qk)(Q-k)) A 2 ^ 2 > 

where = J /ksT is the inverse temperature of the sys- 
tem, and = J2r Q(R)e is the Fourier transform 
of the charge density. 

To proceed further, we shall first replace the Greens 
function mediating the interaction between vortices 
in Eq. J3J with the standard one G^ 2 ' = -V~ 2 
over the square lattice subject to periodic boundary 
conditions 1 ^* 1 ^. This approximation neglects only the 
short-range features of the interaction energy near the 
pinning centers. Secondly, we shall consider the opti- 
mum pinning configuration: V P (R) = — U for points R 
that coincide with the triangular vortex lattice at zero 
temperature, and V p (R) = otherwise. The long-range 
logarithmic interaction among vortices enforces charge 
neutrality with the uniform background charge density 
/, such that: 

J2Q(R) = fL 2 . (4) 
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This means that the system of vortices is incompressible 
at all temperatures. Vacancies and interstitials are there- 
fore impossible at long wavelengths. In the absence of 
extrinsic pins£«i, U = 0, and at low vorticity, / < 1/30, 
the triangular vortex-lattice depins from the underly- 
ing square lattice at a temperature ksTp ^ = 1.5/ J. 
At higher temperatures it "floats", before melting at 
ksTm = J/20. At low vorticity and low tempera- 
ture, the difference in the internal energy between the 
floating and the pinned vortex lattice phases is therefore 
AE = V vx _ p L/, where V vx _ p denotes the average num- 
ber of pinned vortices. The corresponding difference in 
the entropy is AS* = fcsln(2//), which is not extensive. 
The balance in free energy, AE — TAS, yields a (first- 
order) transition between the two phases at a critical 
pinning strength 

U P = N^_ p k B T\n(2/f) (5) 

that vanishes in the thermodynamic limit, L — » 00. Be- 
low we shall confirm this prediction at weak pinning, as 
well as exploring the effect of strong pinning through MC 
simulation of the 2D Coulomb gas ensemble (J2J. 

The following thermodynamic quantities are measured 
in the MC simulations of the Coulomb gas ensemble 
P|l describing a 2D vortex lattice with optimum pins. 
Phase coherence is probed by the superfluid stiffness 
Pinning of the vortex lattice is determined by checking 
for the appearance of Bragg peaks in the vortex density, 
So(k) = \(Qfi}\ 2 - This is quantified by taking the ratio 
between the heights of the first-order and zero-order 
Bragg peaks, which we term the "Bragg ratio" . Intrinsic 
positional correlations among the vortices, on the other 
hand, are measured through the fluctuation contribution 
to the structure function, Si(k) = (Qj:Q_^)~{Q^)(Q_^)- 
Lastly we also measure the six-fold orientational or- 
der parameter characteristic of the triangular vortex 
latticeA-ift 06 = ^ 2 £i,*(exp[6i(0j - 9 k )]). Here Q 3 
is the angle between a fixed direction in space and the 
direction of the bond between the j-th vortex and its 
nearest neighbor. In general, this orientational order 
parameter decays algebraically with system size as L~' 16 
in the thermodynamic limit, L — > 00. The correlation 
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exponent rje is null in the case of a 2D vortex lattice with 
strict long-range order, while r] 6 = 2 if only short-range 
orientational order exist A 

The MC scheme used closely follows that developed 
by Lee and Teiteli&. A single MC move consists of se- 
lecting a lattice point and one of its nearest neighbors at 
random, and adding a unit charge to one of the points 
and subtracting unit charge from the other, thereby keep- 
ing the total charge of the system constant. This alter- 
ation is then either accepted or rejected according to the 
standard Metropolis algorithm. At temperatures lower 
than J/Aks, simulations revealed that the accepted con- 
figurations only possessed charges of Q = +1, the total 
number of charges, N vx , being set by overall charge neu- 
trality (@J. All other configurations, such as those with 
multiple or negative charges, were heavily penalized on 
energetic grounds. This permits a much less computa- 
tionally costly MC move to be employed in this regime, 
which consists of selecting one of these N vx charges at 
random and moving it to an unoccupied neighboring lat- 
tice site. Extensive simulations at low temperatures were 
run using both MC updating methods to confirm that 
they indeed gave the same results, and thereafter the 
second algorithm was used to obtain the bulk of the low- 
temperature results in this paper. The MC estimates are 
based on between 8000 and 12000 measurements, follow- 
ing an equilibration from a random initial configuration 
consisting of 8000 MC sweeps. In all cases this process 
was repeated a number of times, using a different initial 
state, to ensure that the simulations were not becoming 
ergodically trapped. 

Fig. ^ displays temperature profiles of the various 
physical probes that were listed above, for three dif- 
ferent regimes of the strength of the pinning potential. 
The phase-incoherent floating phase lies in between the 
pinned vortex lattice and the vortex liquid as a func- 
tion of temperature in the no-pinning regime shown in 
Fig. 2Ji 8 - 9 . A strange phase- incoherent state that does 
not float, on the other hand, lies in between the conven- 
tional solid and liquid phases in the weak-pinning regime 
(Fig. njj). Table [I] demonstrates that this observation is 
not a spurious size effect. Notice, in particular, how the 
orientational order parameter </>@ and the magnitude of 
the first-order Bragg peak saturate with increasing size, 
and how the phase rigidity remains null throughout. The 
former implies that the "strange" vortex phase exhibits 
strict long-range orientational order, with a correlation 
exponent t/q = Or Fig. [21 displays the resulting phase 
diagram in the U-T plane. The boundary, U P (T), that 
separates the floating phase (I) and the "strange" phase 
(II) extrapolates to zero roughly as L" 1 at fixed temper- 
ature. This was determined from MC simulations at a 
temperature ksT/^nJ = 0.0060 for two different sizes, 
L = 56 and 112, with the addition of the thermody- 
namic limit, L — > oo, under the assumption that U p = 
there. (Note that our MC simulations show metastabil- 
ity in the vicinity of this first-order pinning transition, 



L Ps/J Bragg ratio <f>e 

56 0.001 ± 0.002 0.329 ± 0.061 0.603 ± 0.008 
112 0.000 ± 0.001 0.765 ± 0.003 0.677 ± 0.002 
168 0.001 ± 0.001 0.784 ± 0.003 0.690 ± 0.003 



TABLE I: Finite size study of the Coulomb gas at a tempera- 
ture k B T/2-K,J = 0.0065 and pinning potential U/2nJ = 0.002 
inside the putative hexatic phase. This strength of pinning is 
slightly larger than U p for L = 56, which explains the rela- 
tively small Bragg ratio at this size (recall that U p decreases 
with size). 



which reflects the two-fold orientational degeneracy of 
the floating states.) Consistent with the simple balance 
of free energy (JSJ, we then conclude that the floating 
phase exists only in the absence of extrinsic pins (U = 0) 
in the thermodynamic limit. The last size analysis also 
implies by Eq. (JSJ), that the number of pinned vortices 
in the "strange" phase scales as A vx _ p oc L, which is 
subthermodynamic. Finally, at yet stronger pinning, the 
line U = U m (T) along which macroscopic phase coher- 
ence sets in shows only minor size dependence for the 
three L x L lattices that we simulated. 

A recent theoretical analysis of the 2D Coulomb gas 
(|2*|) in the absence of bulk pinning finds that an inter- 
mediate hexatic phase can indeed exist if rigid transla- 
tions of the vortex-lattice are prohibited 4 . This phase 
contains unbound dislocations that generate appreciable 
fluctuations in the center of mass of the 2D vortex lattice. 
These fluctuations are responsible for both the destruc- 
tion of macroscopic shear rigidity and of macroscopic 
phase coherence at a 2D melting transition. The hex- 
atic phase contrasts with the conventional pinned vortex- 
lattice phase that shows phase coherence. It also differs 
from the vortex liquid phase by the presence of orienta- 
tional order—. We propose to identify the strange inter- 
mediate phase (II) that neither floats nor shows macro- 
scopic phase coherence, but that exhibits strict long- 
range orientational order (see Figs. ^and|21 and TablcQJ, 
with such a hexatic phase. Fig. ^displays the fluctuation 
part of the structure function, Si(fc), in the putative hex- 
atic phase found in the weak-pinning regime. It shows 
six-fold symmetric Bragg peaks of low order, which strik- 
ingly resemble those obtained experimentally in real-life 
hexatic phases 1 ^. Further, Fig. 0] shows a typical config- 
uration of the vortex lattice in the putative hexatic phase 
near the pinning threshold. We may note the presence of 
an unbound dislocation and the small fraction of vortices 
that are pinned. The strict long-range orientational or- 
der (Fig. \3p), the six- fold pattern shown by this intrinsic 
structure function (Fig. |3J) , and the presence of unbound 
dislocations (Fig. 0} are all consistent with the identifica- 
tion of the "strange" vortex lattice with a hexatic phase 
over a commensurate substrate 10 . 

The vortex lattice tends to be fixed to each and every 
commensurate pin at strong U/2ttJ > 0.01 (see Fig. 
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Interstitial vortex/ anti- vortex excitations are the only re- 
maining degrees of freedom in such case. The tempera- 
ture dependence of the phase rigidity @ obtained from 
our MC simulations strongly resembles that of the zero- 
field case (/ = 0) with no extrinsic pins (U — 0)^. In 
particular, Fig. ^ shows how p s decreases smoothly 
from J to at the expected KT transition temperature!!, 

ksTc —^)J- The vortex lattice subsequently frees itself 
from the extrinsic pins at a higher temperature, T p , as 
shown by Figs. Q]and[2| 



III. 3D VORTEX LATTICE WITH OPTIMUM 
COLUMNAR PINS 

We finally apply the above results to the question of 
phase coherence in the vortex-lattice phase of layered su- 
perconductors. In the extreme type- II limit, this system 
can be modeled by an infinite stack of XY-model layers 
with uniform frustration, Eq. JQ), but with an additional 
Josephson coupling, J z — J/7' 2 , between all nearest- 
neighbors across adjacent layers. Also, in this limit 
the magnetic coupling to vortices in adjacent layers can 
be accounted for by weak optimum columnar pinning, 
V P (R), within the "substrate potential" approximation 7 . 
Such a layered XY model can be analyzed through a 
partial duality transformation that is ideally suited to 
the weak-coupling limit. This leads to the following 
partition function that encodes the thermodynamics of 
the coupled system: Z CG = £ { „ 2} (/5 '•/2>y' 2 ) N l n MC\p l ], 
where n z (f, I) is an integer link field on 2D points r 
between adjacent layers I and I + 1—. Here C[p{\ = 
(exp[i ~Y^pi(f)4>{r, I)]) j z= o is the phase auto-correlation 
function of an isolated layer I probed at the dual charge 
that collects onto that layer: pi{r) = n z (r, I — 1) — n z (r, I). 
Also, N[n z ] counts the total number of dual charges, 
n z = ±1. The latter system is dilute in the weak- 
coupling limit reached at large model anisotropy pa- 
rameters, 7' — * 00. It has been shown recently by 
one of us^ that the phase auto-correlator for a pure 
2D vortex lattice that cannot move rigidly has the form 
|C(1,2)| = {p s /J){r Q /r 12 ) mD , with a small 2D correla- 
tion exponent, 772.0 < (287r) _1 . Here tq is of the order of 
the inter-vortex spacing and r±2 denotes the separation 
between the two probes. Since commensurate pinning, 
V P (R), increases phase coherence (see Fig. [Q, the pre- 
vious bound on r\2D continues to hold. Yet the applica- 
tion of the above duality analysis yields a phase rigidity 
across layers equal to! 2 - pjr = (p s /"f' 2 )(r Q l^'a) mD . The 
extremely small bound on r\2D then implies that a de- 
coupling crossover, <C p s /j' 2 , occurs only for astro- 
nomically large anisotropy, fj' 2 > 10 38 . This indicates 
that layer decoupling does not occur in practice in the 
vortex-lattice phase of extremely type-II layered super- 
conductors at temperatures below 2D melting (cf. ref. 
7). 

The above partial duality analysis of course also applies 
directly to the question of optimum columnar pinning in 



the vortex-lattice phase of strongly type-II layered su- 
perconductors. In the strong pinning regime, Fig. \T\c, 
our MC simulation results for a single layer find conclu- 
sive evidence for a standard KT phase transition driven 
by the unbinding of interstitial vortex/anti- vortex pairs. 
This implies that the autocorrelation functions C[p] that 
appear in Zqq are precisely those corresponding to the 
2D XY model in the absence of frustration (/ = 0) and 
extrinsic pinning (U = 0), up to a gauge transformation. 
We thereby conclude that the strongly pinned vortex lat- 
tice at the matching field goes through a superconduct- 
ing/normal transition that is second-order, and that co- 
incides with the universality class of the standard 3D XY 
model. 



IV. DISCUSSION AND CONCLUSIONS 

In summary, we have studied the nature of phase co- 
herence in the 2D vortex lattice at the extreme type-II 
limit with identical commensurate pins through Monte 
Carlo simulation of the corresponding 2D Coulomb gas, 
Eq. J5J). Our most striking result is the identification 
of a strange intermediate phase (II) that is pinned, but 
that shows no macroscopic phase coherence. It lies in 
the midst of a floating vortex lattice, a pinned vortex 
lattice, and a vortex liquid phase (see Fig. |5J|. A recent 
calculation has indicated that an intermediate hexatic 
phase with such characteristics is indeed expected in the 
2D vortex lattice at the extreme type-II limit JTJ if rigid 
translations of the vortex lattice are prohibited 4 . The 
present calculations show that an arbitrarily weak array 
of commensurate pins have this effect in the thermody- 
namic limit [see Eq. (J5J and Fig. ^J. We have proposed, 
in conclusion, that the strange intermediate phase just 
mention is in fact a hexatic vortex lattice. 

On the contrary, theoretical calculations of 2D lattice 
melting in the presence of a weak commensurate sub- 
strate potential do not predict the existence of separate 
liquid and hexatic phases^. It is quite possible, then, 
that the transition that we observe through Monte Carlo 
simulation between a hexatic vortex lattice (II) that is 
pinned and a vortex liquid that is not pinned is in fact a 
sharp cross-over transition. In such case, the phase dia- 
gram shown in Fig. would be topologically equivalent 
to that predicted by theory (see RefiiS, Fig. 4, right-hand 
side). More extensive Monte Carlo simulations may be 
needed to settle this question. 
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FIG. 1: Phase transitions at / = 56 _1 on a 112 x 112 lattice, 
for three values of pinning potential. Thick solid line (black 
dots) = p s /J, thick dotted line (open circles) = <j)Q, thin dot- 
ted line (open diamonds) = Bragg ratio in So(k). Error bars 
are smaller than symbols. All energies are given in units of 
2nJ. 
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FIG. 2: Phase diagram for / = 56 _1 on a 112 x 112 lattice. 
Roman numerals I and II indicate the floating and the hexatic 
phases, respectively. All energies are given in units of 27rJ. 
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FIG. 3: Intrinsic structure function, Si(fc), for the hexatic 
phase (II) at k B T/2ir.J = 0.007 and U/2n.J = 0.002. 



O * • O * O • O O 
"o *0 • • " * 0, • 

# * O , ' • O 0* o» o # o *o 

o o o * o * a • o * o . o ,o\ o* d o * o • 

,0 o o" o" o* o* o r /''o "0 O 0* 
j • O (f cf o o , * o • 

0* # 0* C* Cf \Q/ • o" 0* 
? Cfc # 6 if * , 

O « 0* O O *0 
*J * O g O t , g 0* '0* 

*0 O O, • O , O » O ft 6 * O 

*0 *0 • , , *0 ft *0 £> 

0* 0» • *0 "0 «0 C 

» 0*0*0 *0 *0 * 
o t o*o»o p o * o "o *0 *0 »0 
Q o •jj *jD to o o *n * n * o *q 



FIG. 4: Typical configuration (black dots) for the hexatic 
phase (II) at / = 56 _1 on a 112 x 112 lattice, at a tem- 
perature of kBT/2-Kj = 0.0067 and pinning potential (open 
circles) U/2nJ = 0.0008. Observe the presence of an unbound 
dislocation inside the Burgers circuit shown, and the absence 
of disclinations. Also notice the small number of pinned vor- 
tices (dot in circle) . 



